Point-contact spectroscopy in heavy-fermion superconductors 
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We develop a microscopic model to calculate point-contact spectra between a metallic tip and 
a superconducting heavy-fermion system. We apply our tunneling model to the heavy fermion 
CeCoIns, both in the normal and superconducting state. In point-contact and scanning tunnel- 
ing spectroscopy many heavy-fermion materials, like CeCoins, exhibit an asymmetric differential 
conductance, dl /dV, combined with a strongly suppressed Andreev reflection signal in the supercon- 
ducting state. We argue that both features may be explained in terms of a multichannel tunneling 
model in the presence of localized states near the interface. We find that it is not sufficient to tunnel 
into two itinerant bands of light and heavy electrons to explain the Fano line shape of the differen- 
tial conductance. Localized states in the bulk or near the interface are an essential component for 
quantum interference to occur when an electron tunnels from the metallic tip of the point contact 
into the heavy-fermion system. 



PACS numbers: 74.55. +v, 74.70.Tx, 85. 30. Hi 
I. INTRODUCTION 

Point-contact (PCS) and scanning tunneling spectro- 
scopies (STS) have been widely used to characterize 
the electronic behavior of heavy-electron materials, es- 
pecially the transition into the superconducting state 
and the opening of an excitation gap. Unlike point- 
contact junctions with conventional metals, most heavy- 
fermion materials show asymmetric conductances, which 
have been difficult to explain. In addition to asymmetric 
normal-state conductances, many heavy-fermion super- 
conductors (HFS) are known for strongly suppressed An- 
dreev reflection signalsP^ This is especially true for the 
heavy-fermion superconductor CeCoIns JanlSI Hallmarks 
of superconductivity are phase coherence and Andreev 
reflection ( AR) , which require the existence of a conden- 
sate of Cooper pairs. The AR signature occurs when a 
quasiparticle retro-reflects off a normal-superconducting 
(N/S) interface as a quasihole, while momentum and 
charge get carried across the interface by the Cooper pair. 

The Blonder-Tinkham-Klapwijk (BTK) formulation 
describes the differential conductance in conventional 
N/S junctions remarkably well by invoking a dimension- 
less barrier strength parameter Z, which depends on 
the barrier potential and the mismatch ratio of Fermi 
velocities.^ For HFS this formula predicts that N/HFS 
junctions are in the tunneling limit, i.e., low transmis- 
sion, and AR cannot occur, contrary to experimental ob- 
servations. Attempts to correct this have been ad-hoc 
by postulating boundary conditions at the interface that 
are unaffected by the mass enhancement of the itiner- 
ant heavy electrons.^ Alternatively, a single heavy-band 
tunneling model with an energy-dependent quasiparticle 



lifetime was proposed to explain the strongly reduced AR 
signal, but it lacks to account for the large voltage asym- 
metry of the conductance in the normal stateP^El 

Over time various models for Kondo scattering with- 
out magnetic impurities have been proposed to explain 
point-contact studies ranging from tunneling into two- 
level tunneling systems coupled strongly to the conduc- 
tion electrons^ to localized electron spins at the point 
contact.^ Very recently, Malteseva et alP-^ presented a 
theory for electron cotunneling into a dense Kondo lat- 
tice that can account for a Fano line shape in the con- 
ductance. Also Yanii, :> argued that the PCS spectra for 
CeCoIn 5 are consistent with a two-fluid picture based 
on the Kondo lattice scenario. Common shortcomings 
of all these approaches have been the ad-hoc nature of 
additional parameters to explain the PCS spectra, the 
neglect of localized states at the interface and the effects 
of pair-breaking surfaces in unconventional superconduc- 
tors. 

In this paper, we present a multichannel tunneling 
model for the solution of PCS and STS in heavy- 
fermion materials that is derived from an analysis of the 
PCS measurements of the heavy-fermion superconductor 
CeCoIns. While we focus here on CeCoIns the proposed 
tunneling model has applications to heavy-fermion sys- 
tems in general. For the first time, our multichannel tun- 
neling model quantifies the reduced AR signal and con- 
ductance asymmetry observed in normal-metal/heavy- 
fermion superconductor (N/HFS) junctions without spe- 
cial constraints on tunneling barriers, Fermi velocity mis- 
match, or itinerant band mass renormalization at the in- 
terface. 

The article is organized as follows: In Sec. [IT] we intro- 
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ducc the multichannel tunneling model for the HFS and 



discuss specific limits. In Sec. Ill we present our theo- 
retical results for a PCS junction and compare with sev- 
eral experimental PCS conductance curves varying over 
a wide range of voltage biases, temperature and orienta- 
tions. Finally, we summarize our results in Sec. |IV| 



II. TUNNELING MODEL 

Our theoretical understanding of the normal-state 
properties of heavy-fermion materials is based on either 
the Kondo lattice or periodic Anderson model. The 
Kondo lattice model describes localized electronic mag- 
netic moments at each lattice site coupled weakly to 
an itinerant electron band, while the periodic Anderson 
model describes l ocaliz ed /-electrons hybridizing with 
itinerant electrons! 19 * 20 ! 



A. A microscopic model for heavy fermions 

Solving the Kondo lattice or periodic Anderson model 
is a formidable task for multiorbital materials. Instead 
we model the heavy-fermion materials by two bands of 
itinerant electrons with localized surface states, which 
may be caused by broken /-electron bonds at the surface 
due to the broken translation symmetry, 

a\k,cr ia 

The heavy-fermion Hamiltonian Hhf represents two 
bands of itinerant conduction electrons with band index 
a G {light, heavy} and localized electrons near the sur- 
face with site index i. The operators (£, {ca-ka) create 
(destroy) an itinerant electron with momentum k and 
spin er in band a, while operators // CT (fi a ) create (de- 
stroy) an / electron at site i with spin a. £ a (k) are 
the respective electronic dispersions and E is the energy 
level of the localized / electrons. 

A simple description of a tunneling experiment is com- 
prised of Hamiltonians for the heavy-fermion material, 
the counter electrode, and the transfer or tunneling pro- 
cesses between them: H = Huf + "Hcicctrodc + Wt- The 
counter electrode is given by normal conduction electrons 



^electrode — ^ ^ ^e(^)^fe cr ^fecr; 



(2) 



k,<7 



and the tunneling Hamiltonian describes all possible 
transfers 



— ^2 tk^k'a^kcr^k'a' +tk^-k'a' C a-ka ek ' 
a:k,a\k' a' 



In addition to the standard overlap integrals t a between 
the conduction band in the point contact and itinerant 
heavy-fermion bands there is a finite overlap, ti oc , from 
the point contact to the localized states in the heavy 
fermion. We also account for weak interaction between 
the localized surface electrons and itinerant electrons 
through scattering terms v a (see Eq. 7 below). In gen- 
eral, to get a Fano resonance in the conductance one 
needs interference between different tunneling pathsf^ 
The resulting differential conductance calculated from 
this model Ht will have an asymmetric Fano line shape. 
Our setup is similar to the one for scanning tunneling mi- 
croscope (STM) tunneling through a magnetic atom on 
a metallic surfaceP^EH The difference between an STM 
and a point contact is that while the STM is defined by 
conduction through one or very few quantum channels 
a point contact consists of many conducting quantum 
channels in parallel. Figure [l] shows a cartoon of the 
processes that are active in tunneling between a metal- 
lic point contact and heavy-fermion material. Here we 
will extend the picture of conduction through individ- 
ual quantum channels to a tunneling model to account 
for point contacts on a heavy-fermion material. We con- 
sider strong overlap between electron states in the con- 
tact and the heavy-fermion compound and thus go be- 
yond the strict tunneling limit. We then extract micro- 
scopic model parameters in form of the overlap integrals 
and the energy of the localized state from fits made to the 
asymmetric conductances in the normal state of CeCoIn5 
reported in the experiments by Goll et alW^ and Park et 
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(3) 



itinerant heavy electron slates 



Figure 1: (color online) A cartoon of the different 
tunneling processes from the tip of the point contact to 
the heavy-fermion material (localized and itinerant 
electrons), which are necessary to account for the 
measured asymmetry in point-contact junction 
conductances and reduced Andreev reflection signals. In 
(a) the localized state appears as a broad resonance at 
the interface while in (b) the localized state forms a 
sharp surface state, which acts as a resonant tunneling 
center. 
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B. Tunneling current 



To calculate the tunneling current through a quantum 
channel we employ the s tandard non-equilibrium Green's 
function technique! 26 * 27 ! In Keldysh notation the tunnel- 
ing current per conducting channel is given bj^ 



j(V) = -Trf 3 



tloc ° Gloc,c — t"l oc ° G c Jo 



+ [i*oG ac -i* a oG c 



-i K 



(4) 



The trace (Tr) is a short-hand notation for summation 
over momentum and spin k, a, and the o-product in- 
dicates a folding over common arguments, e.g., U oc o 
Gioc.c = J2k' tioc,kk'Gioc,c(k',k"), and [ } K denotes the 
Keldysh component of the matrix Green's function. No- 
tation for matrices is: a "hat" (x) denotes a Nambu ma- 
trix, while a "check" (x) represents a Keldysh matrix (see 
below). For ease of readability, we suppress the explicit 
dependence of the tunneling elements and the Green's 
function components on momentum and spin. In equa- 
tion |4]) , Gij are Green's function components of the full 
matrix in reservoir space (c=point contact, a= light (I) 
and heavy (h) conduction bands of the HF, and loc is 
the localized state in the HF). It is the components of 
G that straddle the interface (e.g. G c h 7 G c i,G c j oc etc.) 
which enter equation Q. G is to be determined in the 
presence of a voltage bias across the interface. The tun- 
neling Green's function is calculated by summing up sin- 
gle tunneling events in the self-consistent non-crossing 
approximation resulting in the Dyson equation 



G = G°+G oVoG +G oVoG°oVoG 



(5) 



heavy fermion material. CeCoIns is a very pure mate- 
rial with a mean-free path of order of 100 nm. There- 
fore, we assume that the junctions are in the ballistic 
limit where all tunneling events conserve momentum, i.e., 

ta,kk' — t a 5(k k ). 

To further simplify our calculations, we shall assume 
that the itinerant microscopic Green's functions are de- 
scribed by quasiclassical Green's functions near the Fermi 
energy, while it is essential to keep the full energy de- 
pendence of the localized Green's function. In the case 
of non-superconducting electrodes the different unper- 
turbed Green's functions are 



G° = 



/so 



u loc — 



-«A/" c f 3 



e R - e A + 20 c (e)e 



e R -e A + 2(j) HF (e)e 



K 



Vic 



e R 



e A 



e R f 3 - E e A r 3 - E 
1 1 



e R h - E e A r 3 - E 



-A' 



(8) 



(9) 



,(10) 



with e R ' A = e ± ir] and T3 is the third Pauli matrix in 
Nambu space. The distribution functions for the elec- 
trons are <j> c {T) = tanh(s - eV c )/2T, <I>hf{T) = tanh(e- 
cVhf)/2T, and(f>i oc (T) =tanhe/2T. The parameter r\i oc 
gives the fraction of localized states and is a dimension- 
less quantity. The 2x2 Keldysh matrix structure of the 
Green's functions and hopping element above is given as 



e R = 



e K 



1 


1 




1 = 





1 

1 
1 



(11) 

(12) 



or equivalcntly 



G 



(1 - G° o V)- 1 o G° 



(6) 



The transfer matrix V, which is derived from the tun- 
neling Hamiltonian Ht, and sketched in Fig. [T] details 
all possible tunneling processes of electrons between the 
metallic tip, localized and itinerant electrons in the HFS, 
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A: A:' 



loc.kh' 

f* 

L h,kk' 



^loc,kk' ^h,kk 

Vh.kk 

*h,kk' 
ttkk> 



ii,kk' \ 

Vl,kk> 



(7) 
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Here ti OCt h.i-kk' are the wavefunction overlap integrals be- 
tween contact and heavy-fermion material making up 
the tunneling elements between states k and k , while 
Vh,i,kk' ar e surface-induced scattering elements between 
the localized state and the conduction bands h, I in the 



For the Green's function of the localized state, G° oc , we 
assume a single level at e — Eq. To further simplify our 
model, we will assume that only the heavy electrons un- 
dergo a superconducting transition at T — T c , while the 
light electrons remain uncondensed. Similar arguments 
were entertained for the unconventional multiband su- 



perconductor Sr2Ru04,^E2l although one might expect, 
based on general grounds, that both bands go supercon- 
ducting because of interband interactions leading to a 
proximity effect.^ 

The model of a single superconducting band may be 
justifiable for CeCoIns, where no multiband gaps have 
been seen in the PCS data by Park et alPand Goll et alP 
However, thermal conductivity measurements have been 
interpreted in t erms of uncondensed electrons^! or fully 
paired electrons! 33 * 34 ^ In this work, we shall take the view 
point that the light electrons remain uncondensed and 
explore the theoretical consequences for the PCS spectra. 
Thus for temperatures T < T c only the heavy electron 
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Green's function G° becomes 



Gl = N h 



l rf {e)e R + g* arf {e)e' 



, (13) 



with g^ u ff(e) the self-consistently determined surface 
Green's functions. It is worth noting that the sur- 
face Green's functions are calculated in the quasiclas- 
sical approximation. They can fully account for sur- 
face pair-breaking due to the crystallographic orienta- 
tion of the surface, disorder pair-breaking, and realistic 
band struct ureP^HUl Surface pair-breaking can happen 
for anisotropic order parameters depending on the crystal 
orientation, but not for isotropic s-wave order parameters 
for which the standard BTK expressions were derived. 

The g R ' A (e) are the usual (quasiclassical) retarded and 
advanced superconducting Green's functions. Here we 
made the usual quasiclassical approximation of a con- 
stant density of states per unit energy and per spin for 
conduction electrons at the Fermi surface, both in the 
contact of the metal tip (6°) and the heavy fermion 

(G°): ATchM = Yskesz *(e - £c,fc,»(*)) « K,h,l- How " 
ever, we kept the energy dependence of the localized den- 
sity of states. To compactify our notation, we move 
the density of states factors N c ,h,i and i]i oc into the 
tunneling elements by re- writing \Jr\i oc Nc tioc Uoc, 
y/AfcMh,i t>h,i ~> t ht i, y/i]iocNh,i Vh,l -> v h ,i. The new 
tunneling elements th,i are dimcnsionless, while ti oc and 
Vh,i have dimension -^/Energy. 

Now we can solve through matrix inversion for the full 
Green's function G in Eq. With the previous def- 
initions we calculate the Green's function components 
needed in the expression for the tunneling current in 
Eq. Q through a point contact. First, we focus on the 
case of a heavy-fermion material in the normal state. For 
this case, we derive an analytic expression for the differ- 
ential conductance for a single quantum channel, which 
leads to a modified Fano expression for the multiband 
tunneling model 



dl 
dV 



(V) 



^ 1 

V 

h t 



x cosh 



de (<7pT 



E ) 



B 2 



2ir 



T 2 + ( E - E ) 



eV 



2T 



(14) 



In Eq. ( 14 1 D is the transparency, E is the tunneling- 



renormalized position of the localized energy relative 
to the Fermi level, T is the half-width of the reso- 
nance, and qF is the quantum interference parameter 
that controls the resonance shape. The additional pa- 
rameter B is present for multiband models only, when 
tunneling through a resonant localized state couples dif- 
ferently to the HF conduction bands (see below). B 
adds a Lorentzian to the conventional Fano resonance 
which may be absorbed into a complex Fano parameter 



qF — > q = qF + i <Zs with q& = B/T. If one introduces the 
following two angles 9 t and 9 V and writes 



th =t cos 6t, Vh — vcos8 v , 
ti = t sin 9 t , vi — v sin 9 V , 



(15) 
(16) 



then the five phenomenological model parameters intro- 
duced in Eq. ( 14 1 depend on five microscopic parameters 
(i, v, tioc, Eq, 5 = 9 v — 9t) given by the following relations 



V 

r 

B 



4f 2 



= E, 



(1+i 2 ) 2 

2ti oc vt cos S 
1 + t 2 

t2 ioc + v \ l + t2 sin2 S ) 
1 + t 2 

1 {E -Eq)1- t 2 



(17) 

(18) 

(19) 

(20) 
(21) 



The angle 9t is a free parameter that quantifies the rel- 
ative weight of the overlap integrals t^j. 9 t cannot be 
determined from normal-state PCS measurements. The 
transparency of a single channel, T>, depends only on 
^ = \f^i + ^ ne parameter B is non-zero only if 5 ^ 0, 
i.e., 9 t 7^ 9 V . Finally, we see that in order to observe a 
Fano-shaped normal-state conductance, i.e., a finite qF, 
one needs to be in the limit of small to intermediate tun- 
neling coupling, t 2 < 1 or T) < 1, and have a sizable 
renormalization due to tunneling into the localized state, 
Eq 7^ E . 

At this point we like to comment that the differen- 
tial conductance derived in our multichannel tunneling 
model with localized nonmagnetic surface states, shown 
in Eq. (14 1, reduces to the Fano expression by Yang^ 



obtained for a Kondo lattice model with hybridized c-f 
electrons. Hence, we conclude that the observation of a 
Fano line shape in the conductance is not a sufficient con- 
dition for probing bulk / electrons in the normal state. In 
order to associate the Fano peak with bulk / electrons 
additional tests are necessary, e.g., observing the bare 
localized state Eq crossing the Fermi level with temper- 
ature or a magnetic field dependence of the line width. 

The PCS experiments measure a weak Fano-like con- 
ductance superimposed on a large background conduc- 
tance Gq(V). The asymmetry of the Fano shape in the 
conductance with respect to voltage accounts only for 
3-5% of the total conductance. To capture both these 
contributions we need to go beyond the single quantum- 
channel conductance calculated in Eq. ( 14 1. Since a point 



contact forms over a sizable area, S, compared to atomic 
scales (~ kp 1 ), one should expect thousands of channels 
present, each acting as a single quantum channel. The 
conductance should then be written as a sum over the 
contributions of the individual channels 



dl 
dV 



i^channels 



di 

dV 



(22) 
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V(mV) V(mV) V(raV) V(raV) 

Figure 2: (color online) PCS spectra calculated self-consistcntly for a two-band model in the absence of scattering to 
localized states and between bands (i.e. t loc = v = 0). We assumed d-wave superconductivity in one band and 
uncondensed electrons (normal state) in the other band. The black and red lines are for a [100] (0°) and a [110] 
(45°) interface, respectively. We vary the transparency T> from tunneling (left panel) to high transmission (right 
panel). The relative weight of tunneling, 9 t , into either band is varied from top to bottom for given T>. For clarity 
the conductances for different t are shifted vertically by AG. The zero-temperature gap is A =0.6 meV and the 
background conductance is Go = 174 (fcfi) -1 . 



The interface of a point contact is probably not atomi- 
cally smooth, which means that the majority of localized 
states near the surface will be broadened to resonances 
due to destructive scattering in the interface (see panel 
(a) in Fig.^J. In the simplest case, this means that the lo- 
calized state is broadened into a resonance as prescribed 
by E -> Eq + 

ilbroad- If the resonance is broad enough 
the conductance kernel in Eq. ( 14 ) gives only a broad and 



quantum channels becomes 



featureless contribution resulting in the background con- 
ductance. This applies to our analysis if Abroad ~ 20 — 100 
meV is the largest energy scale of the problem. From 
these types of tunneling channels, i.e., those dominated 
by Jbroad, Eq. (14) gives a weakly voltage-dependent 
background conductance Gq(V) ps Go + 5Gq(V). For the 
remaining few channels, the localized state forms a sharp 
surface state with Eq — >• Eq + sharp, where ^ sharp is 
much smaller compared to other energies of the problem. 
To describe these channels we extract a phenomenolog- 
ical parameter Cq = risharpDe 2 /h from the PCS data 
(n s har P is the fraction of channels having a sharp local- 
ized state at the interface). Go gives the proper magni- 
tude of the Fano resonance relative to the background 
conductance. This bimodal distribution of two types of 
tunneling channels can be resolved in the normal state by 
numerical fits to the various data sets being considered 
here. For the tabulated values of Go and Go, see Table I, 
we find that their ratio (Go/Go) is typically between 3- 
5% for the point contacts we analyze. Another important 
point is that both Go and Go are only weakly temperature 
dependent, which indicates that the bimodal distribution 
of channels is a stable feature of each contact realization 
in a PCS experiment. Therefore the phenomenological 
differential tunneling conductance expression for a point 
contact accounting for this type of distribution of single 



dl 

dV 



(V) = Go 



de (q F T + e - E Q ) 2 + B 2 



AT 



T 2 + (e- Eo) 



x cosh 



eV 



2T 



G (V). 



(23) 



The integrand in equation ( 23 1 depends on a set of four 



parameters (Eq, T, qp, B) that can be extracted from PCS 
experiments. These determine four microscopic param- 
eters (tio C (t),v(t), E (t),5(t)) that depend on the trans- 
parency T> through the tunneling parameter t. The fac- 
tors Cq,Gq(V) are determined from the large voltage- 
scale conductance and in principle determine the distri- 
bution of channels (abroad, 1 'sharp)- In what follows we 
will assume the simplest bimodal distribution of only two 
possible values of Abroad > r and 7 'sharp — > 0. So far Q t is 
the only model parameter undetermined by normal-state 
PCS data. 

One should think about the sharp channels as a set of 
single- channel atomic point contacts. If we have a lattice 
mismatch between the tip and the HFS one may expect 
that a two-band model, as the one proposed here, should 
show a distribution of the angle 9 t . It seems natural to as- 
sume that the wavefunction of the tip may have different 
overlaps with each of the itinerant bands in the individual 
atomic point contacts, but still have a uniform transmis- 
sion over the whole contact area. A consequence of this 
will be that different metallic tips, e.g., Au vs. Pt, should 
result in different overlaps with the itinerant bands of the 
HFS. We conclude that a more general formulation of the 
differential point-contact conductance should be written 
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dl\ 
dV 



) = f dd tP (8 t ) 
J PC J 



\dV 



+ G , (24) 



sharp 



where p(0 t ) is the distribution of 9 t and the (dI/dV) s /j arp 
is the numerically obtained conductance. Instead of mod- 
eling the distribution p(0t), we will look at PCS spectra 
for different values of 9 t . 



C. Special limits of the tunneling conductance 



The normal-state differential conductance in Eq. (23) 



has several instructive limits depending on the particular 
choice of microscopic model parameters. We consider the 
following four cases: 

(a) A one-channel tunneling model: Tunneling into one 
band only, i.e., ti oc — and 9 t = 0,7r/2, reduces to 
the standard expression for a single channel con- 
tact with a transparency V = 4i 2 /(l + t 2 ) 2 and 
a featureless differential conductance. The Fano pa- 
rameters qp , I\ and Eq all vanish, thus resulting 
in no Fano resonance. When the single band goes 
superconducting one obtains the standard Andreev 
conductance for HFS. This case corresponds to the 
special limit of 9t = (ti = 0) of the two-channel 
model to be discussed next. 

(b) A two-channel itinerant tunneling model: Tunnel- 
ing into both light and heavy bands but keeping 
Uoc = v — will not generate a Fano resonance in 
the conductance as seen from Eqns. (17)- (21 1 with 
r = qp = Eq = 0. In the normal state and for a 
constant density of states at the Fermi level in light 
and heavy bands, this limit gives a conductance 
which is constant. When the heavy band goes su- 
perconducting, T < T c , the two-band model gives 
an Andreev conductance, which may be reduced in 
signal with respect to the background conductance. 
In Fig. [2] we show conductances for fixed values of 
transparency T>, but with varying relative weight of 
the tunneling elements via the tunneling angle 9 t . 
The key result of these self-consistent calculations 
(for details see Sec. IIIC) is that for high trans- 
parency junctions (D — > 1) and overwhelmingly 
tunneling into paired heavy electrons {9 t < 7r/8) it 
is impossible to differentiate between a rf-wave su- 
perconductor with nodal lines along [100] vs. [110]. 
Note that for 9 t — n/A the re- normalized tunnel- 
ing matrix element for paired heavy and normal 
light electrons is equal. The situation is reversed 
for 9t = 37r/8 when tunneling is predominantly 
between the normal metallic tip and the normal 
light electrons in the HFS, see bottom curves in 
Fig. [2| This generic two-band tunneling model 
demonstrates that in principle PCS data can differ- 
entiate between tunneling preferentially into paired 
heavy electrons versus uncondensed light electrons. 



(c) A two-channel hybridized tunneling model: Simul- 
taneously tunneling into a single itinerant band 
(here we consider a single heavy band, i.e., <; = 
and th 7^ 0) and a localized state i/ oc ^ with 
nonzero scattering between localized and itiner- 
ant heavy electrons, i.e., vi — and v/i 7^ 0. 
This point-contact tunneling setup will generate a 
Fano resonance in the differential conductance with 
qB = B = 0, for details see previous general Eqns. 



( 17 )-( 21 ) . Very recently, Yang® discussed a Kondo 
lattice model of hybridized c-f electrons. He de- 
rived a normal-state Fano conductance similar to 



the one in our microscopic model in Eq. ( 14 1. In or- 



der to fit the Fano resonance in CeCoIn5, he intro- 
duced a voltage-dependent Fano parameter T(V) 
and a large interband term qg > qp, which he at- 
tributed to multiband and correlation effects be- 
yond c-f hybridization of electrons. Assuming a 
constant qg implies a strongly temperature depen- 
dent interband scattering coefficient B ~ T(T), 
which is difficult to reconcile within our model. 
Furthermore, as we have shown above, introduc- 
ing a nonzero is equivalent to introducing asym- 
metric scattering between different itinerant bands 
and localized states, S 7^ 0, invalidating the model 
assumption of only two species of electrons used 
in the standard Kondo lattice model with only c-f 
hybridization. 



(d) A multichannel tunneling model: Simultaneously 
tunneling into localized and itinerant bands is of 
topic interest and will be discussed in the next sec- 
tion. 



III. RESULTS AND DISCUSSION 



We take the following approach for extracting the 
microscopic tunneling elements (t,9 t ,v,ti oc ,5, Eq) from 
point-contact conductance data: We fit the differential 



conductance in Eq. (23) to an experimental dl/dV-curve 



at a given temperature by extracting all model parame- 
ters, i.e., the four Fano parameters Eq, qp,r, B, the rel- 
ative weight of the Fano-like conductance Co, and the 
background conductance Go(V). After extracting the set 
of phenomenological parameters (Eo,qp,T,B) through 
numerical fits, we solve for the set of microscopic pa- 
rameters (Eq, v, ti oc , S) by treating t and 9 t as free fit pa- 
rameters. Consequently, the parameters t or cquivalcntly 
T> = At 2 1 (1+t 2 ) 2 and 9 t can only be determined by study- 
ing the conductance in the superconducting state as they 
effectively drop out from the normal-state conductance, 



see Eq. ( 23 ) 
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Figure 3: (color online) Panel (a): Tunneling fits for several normal-state conductances (o) taken with a Au-tip on a 
[001] surface of CeCoIns for a set of temperatures. The curves are shifted down by 4 (kSl)^ 1 for each temperature 
trace. Panels (b)-(d): Fitted Fano parameters and extracted microscopic model parameters based on dl/dV-curves 
in (a). The dotted lines in (a) are dl/dV-curves calculated from Eq. Q using the microscopic parameters (open 
diamonds 0; closed diamonds are a lower bound for parameters) shown in panels (c) and (d), which were extracted 
from the fitted Fano parameters (solid circles •) shown in panel (b). We fit each model parameter (r and Eq are in 
units of meV) in panel (b) to a temperature dependent function x(T) — x + x 2 ■ (T/45K) 2 and re-calculate dl/dV 
(dashed lines in panel (a)). In panels (c) and (d) we show the dependence of the microscopic parameters (U oc , v) and 
Eq on temperature for transparencies < T> < 0.175. For transparencies T> > 0.175 the fitting procedure fails and no 
solutions are found for Eqns. ( 17 )-< 21 ). Here we set B = 0. 



A. Model parameters in the normal state 

In Figs. [3] and [4| we show results extracted from a set 
of conductances taken at different temperatures using a 
Au-tip on [001]-oriented CeCoIn 5 . For each individual 
dl/dV-curve we get a good fit to a Fano resonance over 
the entire measured voltage window of |V| < 20 mV. We 
note that while the dl/dV-curves are well fitted within 
the chosen voltage window in both Figs. [3] and [4] the 
extrapolated large voltage-scale conductances are very 
different as we assumed a constant background Go in 
Fig. |J while in Fig. g] we modeled the background con- 
ductance as G (V) = G - Gitanh(VyV*), where V* 
is an additional parameter that is always of the order 
of the voltage window and will not be considered any 
further. As seen in these figures, the Fano parameters 
(E (T), q F (T),T(T),B) depend sensitively on the details 
of how the background conductance is modeled. There- 
fore, for deriving meaningful Fano parameters from PCS 
measurements it is very important to measure over volt- 
age biases as large as possible while at the same time 
avoiding heating. 



The model parameters {qp, Eq, T, B, Go, Go) all have a 
temperature dependence that may be fitted with a poly- 
nomial x(T) a £p=r x p ■ (T/T*y, with Pmax = 1 or 
2 and T* is a typical temperature scale for the onset of 
the strongly correlated heavy-fermion state. For CeCoIn5 
we set T*=45 K. When we re-calculate the dl/dV-curves 
with the fitted temperature dependent model parame- 
ters, we find that the goodness of a single fit is sensitive to 
the precise values of qF, Eq, V, and B. These re-calculated 
dl/dV-curves are shown as dashed lines in Figs. [3] and 
|4j whereas the fits with the original parameter sets are 
shown as dotted lines. 

In Figure [5] we show the temperature dependence of 
the model parameters, when PCS data are collected over 
a larger voltage bias window. We discover that (i) the pa- 
rameters depend significantly on the size of the voltage 
bias window over which the fit is performed, although 
they exhibit similar temperature trends; (ii) the inter- 
band scattering parameter B is not uniquely determined 
(not shown). Only the localized level Eq, Fano line half- 
width r, and Fano parameter qp are insensitive to B. 
The transparency parameter Go is weekly dependent on 
B. Since the parameter B has negligible impact on fitting 
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Figure 4: (color online) Tunneling fits for the same data as in Fig. [3J however, assuming a weak voltage dependent 
background conductance SG (V). Panel (a): The conductance curves are shifted down by 4 (kfi) -1 for each 
temperature trace. The dotted lines are dl/dV-curves calculated from Eq. Q using the microscopic parameters 
(solid • and open o circles) shown in panels (c) and (d), which were extracted from the fitted Fano parameters (solid 
circles •) in panel (b). Again, we fit each model parameter in panel (b) to a temperature dependent function 
x(T) = x + xi ■ (T/45K) + x 2 ■ (T/ABK) 2 and re-calculate dl/dV (dashed lines in panel (a)). In panels (c) and (d) 
we show the dependence of the microscopic parameters {U oc , v) and Eq on temperature for transparencies 
< T> < 0.03. For transparencies V > 0.03 the fitting procedure fails over the full range of temperatures and no 
solutions for Eqns. ( IT )- ( 21 ) are found. At low temperatures though, one can find solutions for transparencies up to 
V < 0.047. As in Fig^the parameter B = and A =0.6 meV. 
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Figure 5: (color online) Model parameters vs. temperature from fits to the normal-state conductances taken with a 
Pt-tip on a [001] surface of CeCoIns from Ref. [7] The voltage intervals over which fits were performed was varied 
between \V\ < 20 mV, \V\ < 30 mV, and \V\ < 40 mV, to test for robustness of fit parameters and fit procedure. All 
fits were constrained by setting B — 0. Left panel: Results for a constant background conductance Gq{V) = Go- 
Right panel: Results for a F-dependent background conductance Gq(V) = Go — G\ tanh(F/V*). 



the conductance, we set B = for the remainder of this work. 
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Figure 6: (color online) Tunneling fits to the normal-state conductance data (symbols) at 5 K and 20 K with 
parameters displayed in Fig. [5} For clarity the 20 K data set has been downshifted by 2/kf2. Left panel: Assuming 
Gq(V) = Go, nearly perfect fits are possible in selected voltage windows but not beyond. Right panel: Assuming 
Go(V) = Go — G\ tanh(V/V*), nearly perfect fits are possible in selected voltage windows, as well as outside for 
large voltage windows in the case of the 20 K data set. 
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Figure 7: (color online) Conductance for an s-wave superconductor with model parameters extrapolated from Fig. [3] 
Panel (a): We plot the transparency D{Eq) as a function of the energy relative to the Fermi energy of the localized 
state at the interface. Panels (b)-(d): The conductance is plotted for a range of tunneling angles B t - Panel (e): The 
conductances are compared to the standard BTK conductance vs. transparency V superimposed on the 
normal-state conductance with a Fano line shape. We assumed that 15% of the background conductance is 
superconducting, i.e, t\btk — 0.15 in Eq. (25). Each conductance curve is shifted up or down in steps of 10 (kf2) — 1 
relative to the conductance marked by the arrow. We used Aq =0.6 meV. 



Quite unexpectedly, we find that the magnitudes and 
temperature dependences of the phenomenological model 
parameters depend strongly on the treatment of the con- 
ductance background. The only robust feature that can 
be extracted is a temperature dependent Fano param- 
eter r of order 16 meV that nearly doubles between 5 
K and 30 K. All other parameters depend on the back- 
ground modeling. The extracted temperature depen- 
dence of T(T) suggests the importance of inelastic scat- 
tering at elevated temperatures. Spin-fluctuations seen 
in NMR measurements of CeCoIn 5 and re lated materials 
are likely sources for this behavior! 53 ! 54 ! The renormal- 
ized localized level Eq is either positive (above the Fermi 
level) and nearly independent of temperature (constant 
background Go) or negative (possibly below the Fermi 
level) and decreases further with increasing temperature 
(V-dependent background). We do not observe a cross- 
ing of the localized bare level Eq from positive to negative 
as temperature is increased. This suggests that the local- 



ized states are most likely surface states, since for bulk 
states localized / electron levels should be pushed below 
the Fermi level for temperatures higher than a charac- 
teristic heavy-fermion coherence temperature T* . Above 
T* ~ 45 K the heavy-fermion system CeCoIns exhibits 
well-developed localized / moments. 

To further characterize the origin of the resonance, 
spectroscopic measurements at higher temperatures and 
in magnetic fields may distinguis h betw een localized mag- 
netic Kondo states in a lattice^ESl an d nonmagnetic 
surface states. Currently the picture is not clear. If 
the resonance originates from the Kondo lattice effect, 
one would expect that the resonance disappears above 
the Kondo lattice (coherence) temperature ~ 45 K in 
CeCoIns. However, this characteristic temperature is 
only one quarter of what is expected from the half-width 
of the resonance T ~ fc^T/c ~ 16 meV or Tk ~ 160 
K. Additionally, it has been suggested that a magnetic 
field splits the Kondo resonance due to the Zeeman ef- 
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Figure 8: (color online) Conductance for a d-wave HFS with the same model parameters as in Figs. [| and [7] The 
conductance is plotted for two crystal-to-surface orientations 0° and 45° and for two values of the tunneling angle 9 t . 
The labels nc=non-self consistent, sc=self-consistent, refer to taking surface pair breaking in to account (sc) or not 
(nc). The conductances are compared with the corresponding self-consistent BTK conductance vs. transparency T> 
shown in panel (d). We assume a superconducting fraction i] BTK = 0.15 in Eq. (25). Each conductance curve is 
shifted up or down relative to the conductance marked by the arrow. We used A =0.6 meV. The inset shows PCS 
data from Ref. |H]for comparison. There is a qualitative agreement between the computed PCS in panel (a). 
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Figure 9: (color online) Conductance for an HFS with the same model parameters as in Fig. [4] The conductance is 
plotted for both an s-wave and a d-wave superconductor with two crystal-to-surface orientations 0° and 45° and for 
different values of the tunneling angle t . We used Aq =0.6 meV. The curves have been shifted for clarity. 



feet. However, there is no in dicati on that fields as high 
as 9 T affect the conduct ance ! 1 * 4 ^ 

In Figure[6]we show selected normal-state conductance 
fits for CeCoIn 5 at 5 K and 20 K that were used to ex- 
tract the model parameters shown in Fig. [5] These fits 
highlight the need for conductance measurements over 
voltage biases as large as possible, because fitting a Fano 
resonance over a small voltage window leads to significant 
deviations outside that region and hence quite different 
model parameters. 

In Table|l]we report the temperature dependent model 
parameters for data shown in Figs. [3] and [4j From this 
analysis we find that the only robust Fano parameter 



is r, which measures the half-width of the Fano reso- 
nance, while all other parameters vary from measure- 
ment to measurement and depend on the background fit, 
see Figs. [3j [4] and [5] The temperature behavior of the 
renormalized parameter for the localized level Eq can 
vary from nearly flat to decreasing or increasing with 
increasing temperature depending on the treatment of 
the background conductance. No universal behavior can 
be identified that might relate to the bulk properties of 
CeCoIns, as argued within a two- fluid interpretation of 
the PCS datapSEU How to disentangle surface from bulk 
effects remains a challenge. It should be possible, in prin- 
ciple, to observe the nature of the correlated electronic 
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state and the duality^ of / electrons in CeCoIns with 
tunneling experiments. 



B. Model parameters in the superconducting state 

When a material shows superconductivity one can use 
the nonlinear voltage dependence of its N/S conductance 
to extract further information about the undetermined 
parameters T> (or t) and Q t . In Fig. [7]ja) we plot the 
calculated transparency as a function of Eq for the con- 
ductances fitted in Fig. [3j In panels (b)-(d) we show 
the corresponding conductances calculated at T =0.1 K 
(CeCoIn.5 has a T c ~ 2.3 K) using the extracted temper- 
ature dependent model parameters with an s-wave or- 
der parameter. In panel (e) we show the self-consistent 
BTK N/S-conductance^ superimposed on the normal- 
state conductance (dl /dV)Fano + Go as 



dl 
dV 



dl 

dV Fano 



- Go 



(I-Vbtk)- 



Vbtk dl 

V dV BTK 



(25) 



for the full range of transparencies T>. Here t\btk is the 
fraction of sharp channels in accordance with the model 
for PCS described by equation (23). The self-consistent 



BTK N/S-conductance, (dI/dV)sTK, is computed using 
quasiclassical Green's functions as done in Refs. 1431441 
by accounting explicitly for surface pair-breaking, which 
goes beyond the standard BTK formulation. Therefore 
our self-consistent method enables the determination of 
the spectral properties of the superconducting state by 
incorporating pair-breaking through surface scattering 
for anisotropic order parameters, elastic impurity scat- 
tering, and inelastic scattering off from low-frequency 
bosonic modes. 

In panels (b)-(d) of Fig.[7| we vary the heavy-light elec- 
tron tunneling angle t allowing for competition between 
tunneling into a normal conducting and a superconduct- 
ing band in the heavy-fermion material. As seen, the 
N/S conductance is sensitive to both the transparency 
T>(t) and 9 tl i.e., the relative ratio of tunneling into the 
superconducting vs. the normal-state band. Tunneling 
through a resonant state enhances the effective trans- 



TABLE I: The extracted values of the temperature 
dependent model parameters from data shown in 
Figs. [3] and [4] with constraint B = 0. 



expt. 


Parameter 


Eo 


r 


<7F 


Co 


Go 






[meV] 


[meV] 








Fig. 3 


Xo 


2.01 


13.0 


-2.16 


5.6 


163 




x 2 • (45K)" 2 


1.08 


22.5 


1.59 


2.1 


2.7 


Fig. 4 


x 


1.30 


11.7 


-4.44 


1.66 


166 




xi ■ (45K)- 1 


-14.2 


4.42 


-1.92 


-3.29 


3.81 




x 2 ■ (45K)" 2 


16.4 





5.67 


4.11 


-0.07 



parency of the junction, so that a V(Eq) « 0.17 has an 
N/S conductance similar to a BTK conductance of trans- 
parency V > 0.8. We also see in Figs. [7] and [8] that for 
junctions with "high" transparency T>(Eq) 1 the sub-gap 
conductance is enhanced. Another crucial result of this 
multichannel model is that the conductance enhancement 
due to AR is only ~ 10—15% relative to the normal-state 
conductance and not the conventional 100%. Note that 
the suppression of the AR signal in the HFS comes nat- 
urally about by tunneling into either multiple bands or 
through localized states into one heavy band. Hence the 
reduction of the AR signal can be due to Eq not be- 
ing aligned with the Fermi level and tunneling is slightly 
off- resonant. In turn this leads to an incomplete An- 
dreev reflection as there is a slight parti cle-h ole asymme- 
try with respect to Eq (see G° oc in Eq. (10)). No ad-hoc 
postulates are required for the Fermi velocity mismatch 
between point-contact tip and HFS or special boundary 
conditions at the interface. Instead it is accounted for in 
the tunneling matrix elements of the wavefunction over- 
lap between the metallic tip and the HFS. 



C. Symmetry of the superconducting order 
parameter 

There is quite a body of evidence that the heavy- 
fermion material CeCoIns is an unconventional su- 
per conductor\yrth a d-wave symmetry of the order 
parameter EaUHHIUHlD We follow the strategy outlined 
in Refs. [36 38 43 43] and use the quasiclassical theory 
to compute self-consistently the surface states of a d- 
wave superconductor. The surface Green's function is 
then used to evaluate the conductance taking into ac- 
count surface pair-breaking and hence a reduced order 
parameter at the surface. We emphasize that surface 
pair-breaking of the order parameter is not included in 
the original BTK formulation and hence will lead to dif- 
fering results. The results of our self-consistent calcula- 
tions are shown in Fig. [8] using the same model param- 
eters as in Fig. [7] As seen in panels (a)-(c) we have a 
sensitive dependence on the surface orientation relative 
to the crystal axis orientation.^ If the surface normal is 
aligned with the crystal axis along which the d-wave or- 
der parameter A(T) cos2(</>j, — <f>) has a lobe (</> = 0°), 
then the conductance shows qualitatively the same shape 
as an s-wave superconductor, especially if one takes into 
account the trajectory-average over the order parameter 
A(T) cos2(0j. — <fr). If the surface normal is misaligned 
with the crystal-axis the surface is pair-breaking for a d- 
wave superconductor as the surface scattering connects 
trajectories with different order parameter values, i.e., 
A(T) cos2(0 m A(T) cos2(</w - 4>). For tra- 

jectories where the order parameter c hange s sign there 
is a zero-energy Andreev bound statepEU w hich gives 
rise to a zero-bias conductance peakP^ Therefore the NS- 
conductance in a d-wave superconductor will depend on 
the principal tunneling direction a point contact or STM 
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tip has relative to the crystal orientation. In panel (d) in 
Fig. [8] we show the self-consistent BTK-conductance com- 
puted with Eq. ( p5| ) using the superconducting fraction 
Vbtk = 0.15. For T> < 1 there is a clear dependence on 
the surface-to-crystal orientation. For fully transparent 
contacts on the other hand there is very little difference 
between a 4> = 0° contact and one with (j> = 45°. For the 
dl/dV-curves calculated using the model parameters we 
never reach the fully transparent limit, as for the s-wave 
case in Fig. [7J the resonant enhancement of tunneling 
via a localized state gives for T>(E ) ps 0.17 a sub-gap 
conductance similar to that of the self-consistent BTK 
scenario with T> > 0.8. For smaller T>(Eq) the tunneling 
limit is approached. 

The high transmission case, Z> — > 1, gives good agree- 
ment between the self-consistent BTK conductance and 
experimental data^' at all temperatures below T c irre- 
spective of surface-to-crystal orientation <fi, see Fig. [8] 
Similarly, the conductances calculated with the multi- 
channel tunneling model reproduce to a large extent the 
experimental PCS data, as can be seen in Figs. |H and H] 
In Fig. [9] we show the corresponding conductance curves 
for an s-wave and d-wave superconductor in the pres- 
ence of a voltage-dependent background Go(V), where 
the model parameters are extracted form normal-state 
fits shown in Fig. |4j It is clear from these self-consistent 
calculations that PCS measurements should be able to 
differentiate between tunneling into the nodal vs. the 
antinodal direction of the gap function, irrespective of 
the transparency of the junction. Further, in the Fano 
scenario we find that the amplitude of the conductance in 
the sub-gap region has the proper suppressed magnitude 
compared to the background. The suppression is a direct 
consequence of either competing or interfering tunneling 
channels. 

A key result of these self-consistent calculations is that 
a modified expression for the BTK conductance for an 
HFS point-contact junction, see Eq. (25 1, gives the cor- 
rect description of experiment. However, an unphysical 
parameter T> — > 1, i.e., tunneling in the high transmission 
limit, is required. On the other hand, the multichannel 
tunneling model gives the correct description with phys- 
ically reasonable microscopic parameters, i.e., tunneling 
in the low transmission limit with Pel. 



• The generic two-channel itinerant tunneling model 
demonstrates that in principle PCS data can differ- 
entiate between tunneling preferentially into paired 
heavy electrons versus uncondensed light electrons. 

• A consequence of the self-consistent calculations for 
the two-channel itinerant tunneling model in the 
absence of localized states is that for high transmis- 
sion junctions {D — > 1) and overwhelmingly tun- 
neling into paired heavy electrons (9 t < tt/S) it is 
impossible to differentiate between a c?-wave super- 
conductor with nodal lines along [100] vs. [110]. 

• The only robust Fano parameter, we succeeded to 
extract from several PCS measurements on differ- 
ent samples, is T(T) with T(0) ~ 16 meV. It mea- 
sures the half-width of the Fano resonance. Its 
temperature behavior suggests the presence of sig- 
nificant inelastic scattering, which may be due to 
self-energy effects like scattering off from spin fluc- 
tuations or electrons. 

• Finally, for a multichannel tunneling model the 
zero-bias conductance enhancement due to An- 
dreev reflection is reduced to only ~ 10 — 15% rela- 
tive to the normal-state conductance compared to 
the conventional 100% effect. The origin of this 
suppression can be due to either tunneling into mul- 
tiple competing itinerant bands or through a local- 
ized state into one itinerant band. 

We conclude that it is desirable to have tunneling mea- 
surements ranging from the low to high transmission 
limit. Future measurements at higher magnetic fields, 
higher temperatures, and over wider voltage bias win- 
dows will help to probe quantum interference of electrons 
tunneling between a metallic tip and itinerant heavy- 
fermion bands. In order to identify the origin of the local- 
ized states nanosized PCS tips or STM tips, which are 
in the single-quantum channel limit for tunneling, will 
prove to be critical for resolving the mystery of tunnel- 
ing into heavy-fermion materials. Since our microscopic 
multichannel tunneling model is quite generic, it should 
also be applicable to other heavy-fermion materials. 



IV. CONCLUSIONS 

In summary, we developed a microscopic tunneling 
model for heavy-fermion materials. A narrow spectro- 
scopic feature associated with localized states seen in the 
heavy-fermion material CcCoIns is modeled through a 
multichannel tunneling junction. The asymmetric line 
shape of the differential conductance is understood as a 
Fano resonance for localized states in the vicinity of the 
interface coupled with itinerant heavy electrons in the 
bulk of CeCoIn 5 . 

We list the key results of our analysis: 
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